基于 变 物性 的 格子 - 玻 尔 兹 曼 通 量 求解 器 


摘要 : 本 文 提 出 一 种 考虑 流体 热 物 性 变化 


(中 
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的 影响 ,在 保留 了 传统 格子 - 玻 尔 效 曼 方法 固有 优势 的 同时 ， 殉 服 了 在 网 格 划分 、 边 界 条 件 处 理 等 方面 的 缺陷 。 本 文采 


的 格子 - 玻 尔 兹 曼 通 量 求解 器 。 
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该 求解 器 能 够 捕 提 流体 物性 参数 的 变化 对 流动 和 传 热 


器 模拟 了 同心 环 型 腔 内 由 温 


常 物性 解 低 估 了 换 热 设备 的 传 热 性 能 ， 而 


差 驱动 的 自然 


该 求解 


对 流 , 讨论 了 不 同 温差 条 件 下 的 流动 和 传 热 特性 , 分 析 了 物性 变化 的 影响 。 结果 表明 : 
Ra 数 越 大 ， 偏 差 越 显著 。 
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A Variable Property-based Lattice Boltzmann Flux Solver 


CAO Yu-Hui! 
(1. College of Sciences and Engineering, University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: A variable property-based lattice Boltzmann flux solver is proposed in the present paper. This solver is 


capable of capturing the effects of the variation in fluid properties on flow and heat transfer characteristics. It retains 


inherent advantages of conventional lattice Boltzmann method and overcomes many drawbacks such as the difficulty 


in using non-uniform meshes, the inconvenience of dealing with boundary conditions, etc. In the present paper, the 


natural convection induced by a radial temperature difference in a horizontal concentric annulus is simulated by the 


solver, the flow and heat transfer characteristics obtained under different temperature difference conditions are 


studied, and variable property effects are discussed. It is found that the commonly-concerned constant property 


solution (CPS) underestimates the heat transfer performance of heat exchangers, and the deviation of the CPS from 


the variable property solution becomes more and more notable with the increase of the Rayleigh number Ra. 


Key words: lattice Boltzmann method; lattice Boltzmann flux solver; variable property effects; natural convection 


0 引言 


格子 - 玻 尔 效 曼 方法 (LBM) 基于 介 观 动 理论 ， 
其 中 的 演化 过 程 物理 图 像 清晰 ， 计 算 简 单 ， 容 易 编 


最 近 两 年 ，Shu 等 人 ”提出 一 种 格子 - 玻 尔 兹 


曼 通 量 求解 器 (LBFS), 并 
多 相 流 体 的 流动 与 传 热 


将 其 推广 应 
的 数值 研究 


于 单 相 流 体 、 
073 E LBFS 


中 , 通过 对 格子 - 玻 尔 效 曼 方程 进行 C-E 分 析 得 到 守 


恒 的 微分 方程 ， 然 后 采 月 
进行 离散 化 处 理 , 单元 界 


日 有 限 体 积 方法 对 微分 方程 
上 的 数值 通 量 利用 格子 - 


程 , 计算 是 局 部 的 , 具有 良好 的 并 行 性 和 可 扩展 性 ， 
在 模拟 大 规模 复杂 非 线 性 流动 问题 时 有 具有 独特 优 
H, 已 被 成 功用 于 处 理 各 类 复杂 流动 问题 T, ln 
PERE. mit. SA. SH. Ki 


玻 尔 效 曼 方程 的 局 部 解 进 行 估算 ,有 关 LBFS 的 详细 
介绍 ， 可 以 参考 文献 "。 

变 物性 的 影响 是 过 去 几 十 年 理论 研究 和 数值 
研究 的 热点 问题 之 一 。 已 有 研究 结果 表明 ， 流体 热 
物性 随 温度 的 变化 对 流动 不 稳定 性 有 比较 显著 的 影 
响 " 9， 同时 对 传 热 性 能 的 影响 也 不 容 忽视 。 


然而 , 采用 LBM 研究 这 类 问题 的 工作 ”比较 少见 影响 数值 通 量 的 估算 过 程 。 近 期 ， 本 文 作者 提出 一 
而 且 都 只 能 采用 均匀 网 格 。 种 可 以 同时 考虑 流体 密度 变化 的 算法 (VPLBFS), 在 
本 文 将 LBFS 推广 用 于 变 物性 流体 流动 与 传 热 文献 ”中 做 了 详细 介绍 和 系统 检验 。 


过 程 的 数值 研究 。 关 于 本 文中 的 变 物性 LPS, A — 
Amen. S SBT ama, Sa | TERT BURR SR SOR AES 


系数 、 比 热 容 以 及 体积 力 流体 密度 随 温度 的 变化 。 1.1 宏观 控制 方程 
第 二 ， 可 以 很 方便 地 采用 非 均匀 网 格 。 存储 对 于 变 物性 流体 的 流动 和 传 热 问题 ， 在 忽略 和 儿 性 热 


量 小 ， 只 需 存 储 平衡 分 布 函数 。 第 四 ， 普 适 性 好 ， 耗 散 和 压缩 功 的 条 件 下 ， 控 制 方程 可 以 写成 如 下 形 
方程 中 的 体积 以 及 由 物性 变化 引入 的 附加 项 不 会 式 : 

0,p* V:(pu)-0 (1) 

8,(pu)+V-(puu) - -Vp + V| u(Vu+(Vu)") |-F (2) 

Ta (ð T +u-VT)=V-(KVT) (3) 


: " 1.2 标准 LBE 和 多 尺度 C-E 分 析 
SOLO 完工 , 流 ma EE , 导热 系 
EARKI ALTARS MERE TARRE arabe ahi, 变 物 性 流体 流动 的 标准 让 E 忆 相 


和 比热容 C, Ed APR, F 是 体积 力 ， 只 有 可 以 写成 方程 (4) 的 形式 : 
体积 力 中 出 现 的 流体 密度 是 温度 的 函数 。 


ed =. 
(ieee ee 201... N (4) 
T, 
are pies ; > 2 
其 中， 是 华 标 位 置 ， 忆 是 粒子 运动 的 时 间 步 长 ， ja ypg [ecu lub | (s 
i 7 € ge. 


其 大 小 等 于 格子 的 尺寸 6, e, 是 格子 模型 中 æ 方 


向 上 的 粒子 速度 ，N 是 粒子 速度 的 个 数 ，z, 是 弛 KP, o, 是 加 权 系数 ，c, EER o, Alc, 的 数值 


取决 于 所 选用 的 格子 模型 。 对 于 二 维 问题 ， 最 常用 
的 格子 模型 是 D2qo 模型 ,图 1 给 出 了 D2q9 模型 的 


Ye, for ENSE AYER, fI 是 对 


应 的 平衡 分 布 函数 ， 有 具有 方程 C0 WER: 示意 图 ， 各 个 方向 上 粒子 的 速度 具有 6) 的 形式 : 
0, a=0 
e, =4(cos[(a—1)z/ 2], sin[(a@—1)z/2] )c, a =1,2,3,4 (6) 


2 (cos[(e —5)z / 24 z / A], sin[(& —5)z / 24-z / 4] )e, a=5,6,7,8 


HH, c=6,/6,. M T D2q9 模型 ， 声 速 Os = 0; = 0, = 0, =1/36 o 

c,=c/V3, mA, Aa, = 4/9， FM ae 1) 
p=} f, EXE 

Q,—0,-0,-0,-1/9, a20 a0 

fA SE WLS CD 和 (8): pu-Y e f, = oe f Son 


对 方程 (4) 做 Taylor 展开 ， 得 到 下 面 的 方程 〈9) 
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图 1 D2q9 格子 模型 示意 


Fig. 1 Schematic of the D2Q9 lattice model 


在 多 尺度 C-E 分 析 中 ， 分 布 函 数 f, 以 及 f. o 


时 间 和 空间 的 偏 导 数 按照 以 下 形式 进行 多 尺度 展 
FF; 


fof refi +e fo (10a) 
D a P E (10b) 
Ot On ot, 

V, =, (10c) 


其 中 ，2 是 一 个 小 量 ， 其 大 小 正比 于 努 森 数 。 
将 方程 (10) 代入 方程 (9), 经 过 整理 得 到 下 面 的 
方程 : 


P (p= (4D) 
A (pu) =0 

275 e V (pu) - -Vp* V [uva (Vu) | (12) 
其 中 ， 

u= pe (z, -0.5)6 (13) 


对 于 温度 场 ， 引 入 分 布 函数 8. ， 它 的 平衡 分 布 函 


数 记 作 g“。 在 本 文中 ,对 温度 分 布 函 数 g 采用 了 


fan vr (En V 
2 (ôt 


Chinax ive ERAT 


] 1 P , 
forc k )+0(6t?)=0 (9) 


D2q4 模型 ， 对 8， 满足 的 LBE Xtfr Taylor 展开 和 


CE 分 析 之 后 ， 可 以 得 到 守恒 律 (14)， 宏 观 控制 


方程 (15)， 热 物性 和 弛 丈 时 间 z, 之 间 的 关系 式 


(16). 

r= et=De, ae) 
a” aa 

+V-(ul) =V-(DVT) (15) 

D-,-05. (16) 


1.3 变 物 性 LBFS 的 控制 方程 

从 标准 LBE 出 发 得 到 的 宏观 控制 方程 (11)、(12) 
Al (15) 跟 变 物性 流体 自然 对 流 的 控制 方程 (1)、 
(2) 和 (3) 是 有 差别 的 。 为 了 恢复 方程 (1) - (3)， 
同时 又 能 保留 格子 - 玻 尔 兹 曼 模 型 和 守恒 律 的 简洁 
形式 ， 本 文 将 由 标准 LBE 得 到 的 守恒 律 和 关系 式 直 
接 代入 宏观 控制 方程 〈1) - (3) , 从 而 得 到 下 面 的 
方程 : 


£v (Yr ]-o (17a) 

ot a=0 

9pu EF (17b) 
et 


oT 4 1 D (17c) 
一 +V， "4| 1—-—— je, g"* | | 2 XVTV = E 
ôt PE | 2r Jes: J K 


对 于 二 维 问题 ， 在 变 物性 LBFS A, FRE (17) 
可 以 改写 成 以 下 形式 : 


PW ape (18) 
Ot Ox Oy 

其 中 ， 

W={p, pu, pv, is (19a) 
E={P, I, 11, QY (19b) 


F-(P,n, n, 9,3 (19c) 


Chinai Ade HII 
ChinaX iv& ERATU 


D| 定义 了 一 个 无 量 纲 温差 比 9 ,用 于 衡量 温差 的 大 小 。 
rë =o, 0, r, -wrv2] (19d) 
j 0=(T,-T,)/T, (22) 
P, Fers (20a) 
a 在 本 文 的 三 个 算 例 中 ， 瑞 利 数 Ra 恒 为 105， 普 
Boy e rd (20b) F 7 7 
"A 朗 特 数 Pr EX 0.71, T, 1H73 300K, 6 的 取 值 分 别 
is cal fis (4 ni (20c) 为 0.0101, 0.403 10.8, FFLA, T, 分别 为 303.03K， 
a=0 T, 
n 7 yes s MEX ZEE HE »n 
us a f- ES yn (204) 420.9K 和 540K. PAKKASI, RE u, Ee 
MEE Td WAC, 和 导热 系数 K ， 都 是 温度 的 函数 ， 函 数 表 
= 3 eq y 10 neq (20e) ` 、 
Te 天 í male 达 式 采 用 文献 8 中 的 形式 ， 体 积 力 中 流体 的 密度 
m = Sr 2 LC a (20D 仍 采 用 Boussinesq 假设 。 
yy P ay-ay| Ja 2r, a | 
在 变 物性 LBFS 中 ， 采 用 有 限 体积 方法 求解 方程 SSS | Hl 
(18), 其 中 , 每 个 单元 中 心 的 宏观 物理 量 通过 直接 AI ae | 


求解 离散 化 方程 得 到 ， 相 邻 单元 界面 上 的 数值 通 量 
利用 LBM 的 解 通过 局 部 重 构 进行 估算 。 按 照 有 限 体 
CO — 积 方法 ,对 方程 (18) ERME O, 上 进行 积分 ， 得 — 
v= 到 变 物性 LBFS 需要 求解 的 控制 方程 CD. 一 


T^ dW, 1 
N L=-— Ñ RAS, +F" > R,=(nE+n, F), QD 2， 不 同 温差 条 件 下 稳 态 解 的 比较 。 瑞 利 数 Ra —105 。 


D 


dec MUR EM: 流 线 ， 右 侧 ， 等温 线 。 一 0 0.0101 
Eh, AV ERBE 的 体积 ，AS, 是 第 大 个 交界 -- 4-040. -.- 4-08. 
: M " Fig. 2 Comparison of steady-state solutions obtained under 
面 ( 线 ) 的 面积 KE). R 是 第 个 交界 面 ( 线 ) | : i ' 
different temperature difference conditions. Ra=10°. 
上 的 数值 通 量 。 估 算 R, tt 具体 做 法 ， 请 参考 文献 Left: streamline, Right: isotherm. —— 8, =0.0101; 
-- 0,2040, -:- 0,208. 
“= 8。 得 到 尺 , 之 后 ， 采 用 4 阶 龙 格 - 库 塔 方法 得 到 方 i I . 
T dnd 跟 传统 的 LBM 不 同 , 变 物性 LBFS 可 以 很 方便 地 
7 f (21) 的 数值 解 。 采用 非 均 匀 网 格 进行 数值 求解 ， 而 且 不 必 通 过 分 布 
函数 给 定 边 界 条 件 , 可 以 直接 利用 宏观 物理 量 给 出 。 
2 变 物 性 LBFS 算法 检验 ph don eb oan 
对 于 上 述 条 件 下 空 腔 内 的 自然 对 流 ， 采 用 变 物性 


为 了 检验 变 物性 LBFS 的 可 靠 性 , 在 模拟 同心 环 LBFS 进行 求解 ， 经 网 格 独立 性 检验 ， 发 现 使 用 
型 腔 内 的 自然 对 流 之 前 ， 首 先 重 复 了 文献 “中 的 50x 50 的 非 均匀 网 格 足以 得 出 准确 的 稳 态 解 。 另 
工作 ， 模 拟 了 方 腔 内 的 自然 对 流 。 空 腔 的 上 、 下 壁 外 ， 对 于 边界 条 件 的 处 理 ， 直 接 给 定 四 个 壁面 上 的 
面 是 绝热 条 件 ， 左 、 右 两 侧 分 别 是 温度 恒定 的 高 温 


壁面 和 低温 壁面 , 温度 分 别 记 作 T, AUT, ,方程 (22) 


速度 无 请 移 条 件 ， 左 、 右 侧 的 恒温 条 件 和 上 、 下 辟 化 对 流 场 和 温度 场 都 有 比较 明显 的 影响 ， 随 着 温差 
面 的 绝热 条 件 。 初 始 时 刻 ， 流 体 处 于 静止 状态 ， 流 的 增 大 ， 变 物性 的 影响 逐渐 增强 。 为 了 定量 地 分 析 
体 的 温度 等 于 低温 壁面 温度 。 流体 热 物性 的 变化 对 传 热 的 影响 ， 计 算 了 左 侧 高 温 
图 2 展示 了 不 同 温差 条 件 下 稳 态 解 之 间 的 比较 ， 壁面 上 的 平均 努 塞 尔 数 Nu 。 图 3 给 出 了 不 同 温差 
包括 流 线 和 等 温 线 。 如 图 2 所 示 ， 流 体 热 物性 的 变 条 件 下 得 到 的 Nu 随时 间 的 变化 。 如 图 3 所 示 ， 温 
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差 为 0.0101 时 的 Nu 非常 接近 常 物 性 解 中 的 数值 ，。” 利 数 Ra =50000， 流体 的 物性 参数 ， 如 : BEE u, 
随 着 温差 的 增 大 ， 传 热 速率 逐渐 增 大 。 图 2 和 图 3 | uM - 

比热容 C 和 导热 系数 ， 取 温度 为 T 时 的 值 。 这 
中 展示 的 结果 与 文献 “中 的 结果 吻合 ， 表 明 变 物性 TET * e m 


LBFS 能 有 效 捕捉 到 流体 热 物性 随 温 度 的 变化 对 流 些 定性 的 结果 跟 文献 ”中 一 致 。 在 本 文中 ， 把 常 物 
场 和 温度 场 造成 的 影响 。 性 解 作 为 参考 ， 在 下 文中 将 通过 定量 的 比较 ， 分 析 
3 同心 环 型 腔 内 的 自然 对 流 变 物性 的 影响 。 


本 文采 用 变 物性 LBFS 模拟 了 同心 环 型 腔 内 由 径 向 
温差 驱动 的 流动 和 传 热 过 程 。 同 心 环 的 内 环 半 径 和 : 
外 环 半 径 分 别 记 作 R 和 R 9 TL nj ME OO 
rr=R, IR = 2 .。 内 、 外 壁面 温度 记 作 7 和 7 。 "i. vec M 
跟 第 2 部 分 相同 ， 通 过 不 同 温差 条 件 下 数值 结果 之 i 26 mee ow 
间 的 比较 讨论 流体 热 物性 变化 的 影响 。 了 恒 为 * 
300K, 无 量 纲 温差 比 0 4) FIN 0.0101, 0.403 和 0.8. 
普 朗 特 数 Pr 恒 为 0.71。 初 始 时 刻 流 体 处 于 静止 状 
态 ， 流 体温 度 为 。 速 度 边界 采用 无 滑 移 条 件 。 洛 
径 向 采用 非 均匀 网 格 ， 沿 周 向 采用 了 均匀 网 格 。 
3.1 常 物性 解 

首先 ,采用 变 物 性 LBFS 的 简化 形式 ,忽略 流体 秋 度 、 


Time (S) 
导热 系数 、 比 热 容 随 温度 的 变化 ， 得 到 了 同心 环 型 图 3 左 侧 高 温 壁 面 上 平均 努 塞 尔 数 的 时 间 演 化 。 
腔 内 上 自然 对 流 的 常 物性 解 。 网 格 数 为 180x40 。 图 Fig. 3 The temporal evolution of the average Nusselt number 
4 给 出 了 常 物性 解 中 速度 分 量 和 等 温 线 的 分 布 。 瑞 Re 


u y T 
图 4 常 物性 解 中 的 速度 分 量 和 等 温 线 。 Ra-50000 


Fig. 4 Distributions of velocity components and isotherms in constant property solution. Ra = 50000 


3.2 变 物 性 的 影响 为 了 定量 地 分 析 物 性 变化 对 传 热 速率 的 影响 ， 
在 这 部 分 , 采用 变 物 性 LBFS 模拟 了 三 个 不 同 0 条 件 


我 们 计算 了 外 圆 环 上 的 传 热 系 数 K ,按照 文献 四， 
下 同心 环 型 腔 内 的 自然 对 流 。 流 体 的 秋 度 ， 导 热 系 M i 


A ERA A A TE BY PRA, KARRIKA i, DLE XN (23) 的 形式 

8) 中 的 形式 。 Ra 和 Pr 分 别 为 50000 和 0. 71. 
图 5 展示 了 不 同 0 条 件 下 数值 解 之 间 的 比较 , 包 de oun rr -ln(rr) Lao 

括 速 度 分 量 和 等 温 线 的 分 布 . 跟 方 腔 中 的 情况 类 似 ， " 2n(rr-D)^ ôr 

热 物 性 的 变化 对 速度 场 和 温度 场 都 有 比较 明显 的 影 

响 。 而 且 ， 随 着 温差 比 的 增 大 ， 物 性 变化 的 影响 逐 

渐 增 强 。 


(23) 
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u y T 
图 5 变 物 性 解 中 速度 分 量 、 等 温 线 之 间 的 比较 。 — 6,200101, sss 6-040, -.-0 =08。 
Fig. 5 Comparison of velocity components and isotherms in variable property solutions. —— 6, 20.0101, === 0, 20403 , 
eve d x8. 
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Time (s) 


图 6 外 圆 环 上 传 热 系 数 随时 间 的 变化 。 


Fig. 6 The temporal evolution of heat transfer coefficient 


calculated on the outer wall. 


图 6 给 出 了 不 同 温差 条 件 下 k,, 随时 间 的 变化 。 


eqo 


248 =00101 时 ， 物 性 变化 相对 比较 小 ， 此 时 的 变 物 


性 解 与 常 物性 解 之 间 的 相对 误差 很 小 ，K,,, 略 大 于 


常 物性 解 中 Ko 的 值 。 随 着 温差 比 的 增 大 ， 物 性 变 


eqo 


化 的 影响 逐渐 增强 , 当 9 = 08H], ZWEE K 比 


eqo 


相同 Re 和 Pr 条 件 下 常 物性 解 中 的 Kk,, 大 15%。 这 个 


eqo 


定量 的 结果 与 空 腔 内 的 结果 比较 接近 。 这 表明 : 一 
方面 ， 变 物性 的 影响 对 设备 几何 结构 的 变化 并 不 敏 
感 ， 男 一 方面 ， 在 数值 计算 中 ， 通 过 忽略 流体 热 物 


性 的 变化 而 得 到 的 常 物性 解 低估 了 对 流 换 热 系统 的 
传 热 性 能 ， 而 且 温 差 越 大 ， 由 热 物 性 变化 引起 的 偏 


10000 100000 
Ra 


图 7 不 同 瑞 利 数 条 件 下 外 圆 环 上 传 热 系 数 的 比较 。 


ls 


Fig. 7 The heat transfer coefficient obtained under different 


Rayleigh number conditions. 


为 了 进一步 分 析 变 物性 的 影响 ， 本 文采 用 变 物 性 


LBFS 模拟 了 不 同 


Ra 条 件 下 同心 环 型 腔 内 由 温差 驱 


动 的 自然 对 流 ， 计 算 了 稳 态 阶段 外 圆 环 上 的 传 热 系 


Be. ATH 


eqo 


Hf ^ Ra 下 变 物性 解 中 的 K 


eqo ? 


由 图 可 见 ， 随 着 Ra KNK, K, 逐渐 增 大 ， 流体 热 


物性 的 变化 对 


同 温差 条 件 下 ， 
越 大 。 


eqo 


J 影响 逐渐 增强 。 换 言 之 ， 在 相 


Ra 越 大 ， 变 物性 对 传 热 的 影响 也 
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45 论 

本 文 提 出 一 种 基于 流体 物性 变化 的 格子 - 玻 尔 
兹 曼 通 量 求解 器 〈( 变 物性 LBFS )， 该 求解 器 保留 了 
传统 LBM 的 优势 , 同时 克服 了 传统 LBM 在 网 格 划 
分 ， 边 界 条 件 的 处 理 等 方面 的 缺陷 。 有 关 同 心 环 型 
腔 内 由 温差 驱动 的 自然 对 流 的 数值 研究 结果 表明 : 
通常 关注 的 常 物性 解 低估 了 换 热 设备 的 传 热 速 率 ， 
尤其 是 在 温差 较 大 ， 瑞 利 数 也 比较 大 的 情况 下 。 
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